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Abstract: It has recently been pointed out that the usual procedure employed in order 
to obtain the transverse-traceless (TT) part of metric perturbations in lattice simulations 
was inconsistent with the fact that those fields live in the lattice and not in the continuum. 
It was claimed that this could lead to a larger amplitude and a wrong shape for the 
gravitational wave (GW) spectra obtained in numerical simulations of (p)reheating. In 
order to address this issue, we have defined a consistent prescription in the lattice for 
extracting the TT part of the metric perturbations. We demonstrate explicitly that the 
GW spectra obtained with the old continuum-based TT projection only differ marginally 
in amplitude and shape with respect to the new lattice-based ones. We conclude that 
one can therefore trust the predictions appearing in the literature on the spectra of GW 
produced during (p)reheating and similar scenarios simulated on a lattice. 
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1 Introduction 

Gravitational waves (GW) are expected to be produced copiously in the early universe in 
processes like (p)reheating after inflation [1-7], phase transitions [8-11], during the turbu- 
lent motion of plasmas [12] and from the self-ordering dynamics of Goldstone fields [13, 14]. 
As opposed to the GW background generated during inflation from the quantum fluctua- 
tions of the metric [15], these post-inflationary processes correspond to classical mechanisms 
of GW production, due to the motion of large overdensities. 

Such backgrounds of GW could open a new window into the early universe by providing 
precious information about the dynamics responsible for its production, much before pri- 
mordial nucleosynthesis. The very violent dynamics of the fields sourcing the GW cannot 
be described in linear perturbation theory, and usually takes place very far from equi- 
librium. This is the reason why the study of the GW production in the early universe is 
usually done with the help of lattice simulations. It is therefore of great importance to have 
a precise handle on those simulations to be sure that the predictions made on the amplitude 
and shape of the GW spectra is accurate enough for the future GW observatories [16] to 
detect and constrain these backgrounds. 

In homogeneous and isotropic background spaces, the six (independent) physical de- 
grees of freedom of the metric split into irreducible representations of rotations. There are 
two scalar, two vector and two tensor perturbations. The two tensor components corre- 
spond to the two polarizations of the GW, the transverse and traceless (TT) degrees of 
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freedom {d.o.f.) of the metric perturbations. The flat Friedman-Robertson- Walker (FRW) 
line element with metric perturbations in the synchronous gauge, can be written as 

ds^ = -dt^ + a^{t){5ij + hij{yL,t)^dx'dx^ , (1.1) 

The equations of motion of TT d.o.f. of hij are 

d^d^'hj^ = 16ttGUJ^ , (1.2) 

with n^-"'" the transverse-traceless part of the full anisotropic-stress tensor Ilj-,-. 

The transverse-traceless tensor tensor Uj^^ is obtained by applying a projector Aj^^m 
in momentum space 

n,7(k, t) = Aij,iM fiUK t) (1.3) 

where 

AijM^) = PiKk)Pim(k) - ^Pij(k)PiM , (1.4) 
-Pij(k) = dij kikj , (1-5) 

and k = k/|k|. 

Because the projector is non-local in space, applying it to the source 3-t every time 
step is computationally expensive. In practice, it is therefore more convenient and still 
mathematically equivalent (see [4] for details) to consider a tensor hij which satisfies the 
equation of motion 

d^,d''hi, = 167TGUij, (1.6) 

and apply the projector 

hj^{k, t) = AijV^(k) himiK t) (1.7) 

only when calculating the output. Fourier transforming back hjj^{]<., t) to coordinate space, 
one finds that the metric perturbation 

/iT.T(x,t) = j j^e-'^^^Aij^iMhUKt) (1.8) 

verifies the required conditions 

= hj^, (Symmetry) (1.9) 

^/jTT^Q^ (Ti-acelessness) (1.10) 

i 

^ Vi/i,^"^ = 0, (Transversality) (1.11) 

i 

at all X, t, necessary for identifying hj^{'x.,t) with the gravitational wave d.o.f. It is also 
easy to check that the projector (1.4) is maximal in the sense that it leaves any tensor cif^ 
that satisfies the conditions (1.9)-(1.11) unchanged, 

J^ijM'^Jm = <^lj^ ■ (Maximality) (1.12) 
This guarantees that we capture the whole transverse-traceless component. 
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The above procedure to obtain the TT d.o.f. is well-defined in the continuum. How- 
ever, it has recently been pointed out in Ref. [17] that on a lattice one needs to pay more 
attention to the definition of the projector (1.4), since one can define the momentum k in 
many different ways. In particular, Ref. [17] claims that if one applies the wrong projector, 
a significant leak of scalar modes into the tensor modes (GW) might occur, significantly 
modifying the amplitude of the GW spectrum. In the context of (p)reheating, all lattice 
simulations carried out in recent years by the different groups [3-7], see also Refs. [18-20], 
filtered the TT metric d.o.f. with the projector defined above. Therefore, whether such 
projector is or not appropriate for the lattice, challenges the validity of the GW spectra 
shown in the literature. 

In this paper we investigate this issue in detail. We will review first, in Section 2, some 
ideas about the discretization aspects in a lattice, and then in Section 3 we will present 
our method for obtaining a TT-projector consistent with the symmetries of the lattice. 
In Section 4 we will compare analytically and numerically several GW spectra obtained 
with the continuum-based projection and with different lattice-based projections. Finally 
in Section 5 we will summarize and conclude. 



2 Lattice discretization 
2.1 The lattice derivative 

When one simulates the dynamics of non-equilibrium fields like in (p) reheating, the field 
equations are discretized on a lattice. We consider a lattice with A^'^ points (representing 
spatial comoving coordinates) labeled as n = (ni,n2,n3), with nj = 0,1,. ..,iV — 1. A 
function /(x) in space is represented by a lattice function /(n) which has the same value 
as /(x) at X = n5x, with Sx = L/N the lattice spacing, and L the length of the lattice. 

To discretize the equations of motion, one has to replace the continuum derivative with 
a lattice expression that has the same continuum limit. A simple and symmetric definition 
of a lattice derivative is the neutral one, 

|V»/l(n) = + 

where i is the unit vector in direction i. This has the drawback that it is insensitive to 
variations at the smallest scale of one lattice spacing. In this sense, a definition involving 
the nearest neighbors may be preferable. A common way to do this, is to define the forward 
and backward derivatives 

[vf/i(„) = ^n-^^l^n-) ^ 

but these definitions lack the symmetry of Eq. (2.1). This issue can be solved by defining 
the derivative at half-way between the lattice sites, as 

[V./](n + £/2) = ^(" + ;)-^^") . (2.3) 

ox 

To improve accuracy, one can also consider lattice derivatives which involve more points, 
but in practice the definitions have a symmetry either around a lattice site as Eq. (2.1) or 
half-way between lattice sites as (2.3). 
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In order to extract the transverse-traceless component of hij, one needs to apply the 
lattice version of the projector (1.4). Since the projector is defined in Fourier space, on the 
lattice one has to use the discrete Fourier transform (DFT) /(n), defined as 

= W^ "^"^"^ ' = ^ e+¥-/(„) , (2.4) 

n n 

where the index ii = (ni, n2, na) labels the reciprocal lattice, with = — y + 1, — ^ + 2, ... 
— 1, 0, 1, Imposing periodic boundary conditions in coordinate space, i.e. f{n+iN) = 
/(n), there will be necessarily a minimum infrared (IR) momentum kj^ = ^ in the Fourier 
space, such that n will be representing the continuum momentum k = (n-i, n2, na) A:ir. 
Consequently, there will also be a maximum ultraviolet (UV) momentum k\jy = ^^IR P^r 
dimension, and periodic boundary conditions like /(n + iN) = /(ii). 

In the calculation of GW production, there are four separate places where one needs to 
take care of the discretization details of derivatives: in the equations of the fields sourcing 
the GW, in the two sides of Eq. (1.6), and in the transverse-traceless projection (1.7) itself. 
Ideally one should use a consistent choice of a lattice derivative everywhere, but in some 
cases there are restrictions that make this very difficult, for example when dealing with 
gauge fields [7]. It is enough in any case to have a consistent choice among Eqs. (1.6) 
and the equations for the GW sources, whilst the discretization details in the transverse- 
traceless projection can be considered separately. In this paper we focus precisely on such 
details of the transverse-traceless projection in a lattice. 

2.2 The transverse-traceless component of the metric 

With the transverse-traceless projection, we want to obtain a tensor hj^^ that satisfies 
the three conditions (1.9)-(1.11) on the lattice, with respect to the appropriate lattice 
derivative. In the literature^, this projection has been done by taking Ajj^/m(ii) as in 
Eq. (1.4) evaluated at k = iifciR, 

/^^(n) = ]^ E e-'^""A.,7™(n) himin). (2.5) 

n 

However, as highlighted in Ref. [17], the resulting quantity is not transverse with respect 
to any of the usual lattice derivative operators Vi, i.e. VihJ^^ ^ 0. 

In particular, due to this lack of transversality, Ref. [17] claims that a significant leak 
of scalar modes into the tensor modes might occur in such a way that the amplitude of the 
GW spectrum extracted with the above continuum-based projector (2.5) could be several 
orders of magnitude higher than it should be. 

We therefore need a general and consistent procedure in order to define a TT-projection 
in the lattice, i.e. a projector Ajj^/m that restores the transversality with respect to a given 
lattice derivative. Only then we will be able to quantify the potential distortion of the 
GW spectra with respect to the results obtained with the continuum-based projector. In 
order to construct the lattice equivalent of Eq. (1.4), we need a lattice momentum k. Such 
momentum will depend, of course, upon the choice of the lattice derivative with respect to 
which the transversality condition will be attained. 

""^Some papers, for instance Ref. [3] and Ref. [5], considered the projection at the level of the source, 
whereas others like Ref. [4], [6] and [7], considered the projection at the level of the metric perturbations, 
as we are discussing here. 
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2.3 The lattice momentum 

The lattice momentum is given by the Fourier transform of the lattice derivative Vj. To 
keep the discussion general, we do not assume for the moment a specific form for the 
derivative, but simply assume that it is given by a linear operator in the space of lattice 
functions. Therefore, the value of the derivative [Vj/] is a linear combination of the field 
values at different lattice sites, 

[V,/](n)= J]A(n,m)/(m), (2.6) 

m 

where Di{n,m) is a real- valued function of two variables on the lattice. For example, for 
the neutral derivative (2.1), we have 

A(n,m) = — . (2.7) 

Because we want the derivative to be translation invariant, Di{n,in) is only a function of 
the difference n — m, i.e. Di{n, m) = Di{n — m), and we can write 

[Vif] (n) = ^ A(n - m)/(m) = A(m')/(n - m') . (2.8) 

m m' 

For the neutral derivative (2.1), we have 

A°(m') = (2.9) 
For the nearest-neighbor derivative (2.3), m' is half-integer, and one finds 

A(mO = ^"^ '-"/^ (2.10) 

Ox 

More generally, any odd function with compact support will give a meaningful definition 
of a lattice derivative. 

The Fourier transform of the derivative Vj/ is 

vTfin) = e'^"-"[V,/](n) = ^ e^-- ^ D,{n - m)/(m) 

n n m 

= ^ e"^-"' A(n') e'^"-'"/(m) = -.keff(fi)/(fi) , (2.11) 

n' m 

where the effective momentum kcfr(n) is given by 

keff(ri)=i^e'i^-nA(n). (2.12) 

n 

Conversely, any function kefr(n) that has the correct leading behaviour in the Taylor ex- 
pansion of the IR limit |n| <^ N, i.e. kefr(n) « ii/ciR, defines a lattice derivative through 
the inverse Fourier transform. 

For example, the neutral derivative (2.1) gives 

te°... = ""P;^""'^' . (2.13) 
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The forward/backward derivatives (2.2) give 



,± _ o ±i7rn,/jv sin(7rrai/A^) _ s\.u{2TTfii / N) _^ . 1 - cos{2mi/N) 



and the symmetric nearest-neighbor derivative (2.3) gives 

^e.. = 2?^^^(^. (2.15) 

In general, if the lattice derivative is anti-symmetric, i.e. Z)j(— n) = — I?j(n), then the 
lattice momentum keg is real. This is the case also for the derivative used in Ref. [17]. 

3 The transverse-traceless (TT) projection on the lattice 

In this Section we will define a lattice projection operator that satisfies the condition (1.9)- 
(1.12), and which therefore gives the TT d.o.f. of metric perturbations living on a lattice. 
Since the transversality notion on a lattice is associated to the choice of a lattice derivative 
Vj, we will introduce a projector that will guarantee tracelessness and transversality with 
respect to any Vj chosen. 

3.1 A real TT-projector 

Let us start with the simpler case of a real momentum, for example k^g in Eq. (2.13). In 
this case, the projector can be defined in the same way as in continuum. 
In analogy with Eqs. (1.4) and (1.5), we define 



and 



AO.,„(fi) ^ /^0(fi)F,o„(ri) - ^4(ri)Pl(fi) . (3.2) 



Using the properties 

fcOff,^0(n) = 0, /^0(fi)/^0(n) = /^0(n), Pf^{^) = P^,{t,), (3.3) 

it is then straightforward to prove that hj^ {iv) = A^j i^{h)himii>-) satisfies the required 
conditions (1.9)-(1.12): 

Symmetry: 

AO^_,^(fi) = AO.„,(ri) ^ /iJT(fi) = ^TT(~)_ ^ /jTT(n) = /iTT(n). (3.4) 



Tracelessness: 



AL„^(n) = ^ hnn) = ^ /ir(n) = 0, Vn (3.5) 
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Trans versality: 

A:Off,,A°. ,^(fi) = fe°ff/S^(n) = V>;^.T(„) = 0, V n (3.6) 

It is also easy to see that the resulting tensor hj^ is real in coordinate space. If hij{n) 
is real, its Fourier transform satisfies h*j{n) = /ijj(— n), and then because A^j ^^{h) is real 
and even, 

hj;^*{n) = Ali^{h)h*^{h) = AO.,^(-n)/i,,(-n) = hj;'{-h) h[;^*{n) = hj^{i^). 

(3.7) 

Finally, to prove the maximality (1-12) of the operator, we assume a tensor af^ that 
satisfies the conditions (3.4)-(3.6), and note that then it also satisfies Pfja^k = ^ijO^Tjk- 
Therefore we have 

aO „TT _ / s: r ^pOx 1 ^TT _ „TT lpO„TT_^TT (n o\ 

^Hj,lm"'lm — [OiiOjm — ^J^^ijOlm I "Jm " "ij ~ 2 " ~ ' ^ ' 

Of course, all properties just discussed apply, not only to k^^ in Eq. (2.13), but to any 
lattice momentum keff as long it is real. The case of a derivative with an associated real 
lattice momentum, is therefore a simple generalization of the continuum case. 

3.2 A complex TT-projector 

In the more general case, the lattice momentum is complex. For example, this is the case 
with the forward/backward derivatives (2.2) and the associated momenta k^. Thus we 
will be forced to take a projector that is also complex. 

Thus, we look for a projector Pij that satisfies fceff,i-Pij(k) = 0. In order to do this, we 
define 

P^Jii^)=S^,-^-^f^, (3.9) 

with |kegp = kes^i*kcs^i. This projector is complex and satisfies 

1) keS,iPij(k) = 

3) keSJP^J{k)^0 

4) /^*(fi)=P,,(fi) 

7) Pi,{h)Pji{h) = Pain) 

In words, this projector is Hermitian, symmetric under Parity transformations n f-)- — n, 
transverse to k but not to k*, and idempotent (P^ = P) but with no inverse (^ P^^) and 
non-idempotent modulus (PP* ^ P). Demanding property 1) we arrived at the form (3.9), 
and then properties 2) — 8) simply followed from such form. 

If Ajj^/m(n) was built as in the real case (3.2), then property 3) would prevent hj^^{n) = 
^ij,imii^)him{i^) from being traceless. We are thus forced to redefine also A^ym in order 
to guarantee the desired TT properties. Moreover, since Ptj is now complex, so is Ajj 
Therefore we must also ensure that hj^{n) = DFT{Aijim{'n)hi,n{ii)} is real. From the 



2) A;:ff,/^,,(k)/0 

^)k:s,P^J{k)=0 



6) Pij(-h) = Pji{h) 
8) Pij{n)PiJh) ^ Pu{n 



(3.10) 
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properties of the Fourier transform and demanding hj^^ (n) = hj^^{n), the latter condition 
can be achieved if and only if Ajjym(n) satisfies 

A^,^(n) = A,,-,„(-n) (3.11) 

This condition suggests how to build the new projector. We can define 

Ai.Mi^) = Pu{h)P*^{n) - ^P,,(n)Pr^(n) , (3.12) 

which verifies 

Aij,im(n) = Ay-i^(-fi) = Aji,mi{ii-) = A*^Lji{ii) = Mm,ij{n) , (3.13) 

From here it is easy to prove that hj^^{n) = DFTjA^. j,^(n)/i;m(n)} is traceless and trans- 
verse, as well as real. For completeness, let us show explicitly how we obtain these condi- 
tions using (3.12) and its properties: 

Tracelessness (1.10): 

hrih) = Pu{iv)PUiv)hUrv) - ]^Pu{iv)PUn)him{i^) 

= Pmi{r^)hira{ii) " PL{ii)him{i^) = (3.14) 

Transversality (1.11): 

keS,ih]:^{-n) = kcs,iPii{n)P*^{n)him{n) - ^kcs,iPij{n)PI'^{n)hi„,{h) 

= 0-0 = (3.15) 

Reality: 

h]^\n) = ^e+^^^-^-""A^,^(n)/iL(n) = ^ e+^^'H5-n^*.^^^(n)/,,^(_n) 

h h 

= ^e-''^^'-^'^A*^^U-i^)hUn) = ^ e'^'^^^'^-^^A,,- ;^(n)/i,„(fi) = hj;^{n) 

ii n 

(3.16) 

However, we find that the resulting tensor hj^^ is not symmetric , i.e. /i^-^(n) 7^ hj^{ii). 
Similarly, the following properties, which distinguish between the first and the second index 
of hj^{n), are also verified 

1) /cefr,^/i^-^(n) = , 2) kl^^^J^'iii) / , 

3) A;eff,/.TT(„) ^ , 4) fc:ff,/.^.T(n) = . ^'-^'^ 

All these asymmetry aspects are simply a consequence of the properties l)-4) of Pjj, listed 
above, which reflects the fact that Pij is not symmetric but rather Hermitian. 

Related to these issues we also encounter a subtle aspect about the maximality con- 
dition (1.12). We find that Ajj^im(n)j4j^(n) = A^^ {n) only holds for those transverse- 
traceless symmetric rank-2 tensors which are transverse, not only with respect the 
lattice derivative Vj (with lattice momenta k^.^ used to build Aij^i^), but also with respect 
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the conjugate derivative defined through the lattice momentum k*g. For example, if 
one builds the projector (3.12) with the lattice momenta associated to forward derivatives 
V+, then Aij, 

im is only maximal with respect those tensors which are transverse both with 
respect to forward and backward derivatives, i.e. Vf Aij = S/~ Ay = 0. 

In the IR limit |n| <^ N, both /cefr,j and k^g ^ approach the same momentum, nfeiR, 
and thus the full maximality condition and the symmetry under the exchange i o j, 
are recovered. Thus for arbitrarily big lattices these caveats should not be relevant. In 
reality, we are of course limited by computer memory and the lattice sizes we can typically 
consider have no more than = 128, 256 or 512 points per dimension, depending upon the 
field content. Nevertheless, despite these two caveats about the maximality and the even- 
symmetry, the projector defined by eqs. (3.9), (3.12), is one which generically guarantees 
reality, transversality and tracelessness on a lattice, and recovers maximality and even- 
symmetry in the IR limit. Thus, any GW spectra obtained by this method should be 
reliable at least in the IR region of the Fourier space. 

3.3 General projector 

To understand the difficulties faced in the complex case, let us note that the conditions of 
symmetry (1.9), tracelessness (1.10) and transversality (1.11) on a lattice involve comparing 
and adding together different components of the tensor hj^^ and its derivatives. For this 
to be meaningful, these components should be arranged in a symmetric way on the lattice. 
This is not an issue for derivatives defined on lattice sites, such as the neutral derivative 
of Eq. (2.1), which is why the TT projector defined with real lattice momenta verifies 
nicely all of the required conditions. For those derivatives defined halfway between two 
lattice sites, with a complex lattice momentum, it is however something to take care of. 
To illustrate this, let us consider the forward derivative V^. The transversality condition 
becomes 



which is not symmetric under parity, because it involves neighboring points only in the 
positive directions. 

Thus, instead of separate asymmetric Vj and V* derivatives with complex lattice 
momenta, such as the forward and backward derivatives (2.2), we should use the symmetric 
version (2.3). In order to make the lattice transversality condition symmetric under parity, 
we should then define the tensor hj^ not on the lattice site n, but at the point n+£/2+j/2, 
which corresponds to the center of a plaquette spanned by unit vectors i and j starting 
at point n. More precisely, off-diagonal components are defined at plaquettes, but the 
diagonal components, for which i = j, live on lattice sites. To avoid confusion, we denote 
the tensor defined in this way by h^"'"(n + i/2 + j/2). 

With this definition, the symmetry of the tensor. 



is a meaningful concept because the two sides of the equation are defined on the same 
plaquette n + t/2 + j/2. The same is true for the trace hjj(n) because all the terms in the 
sum are now defined at the same location. 




(3.18) 



hj> (n + i/2 + j/2) = hlj' (n + i/2 + j/2) 



(3.19) 




(3.20) 
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and similarly in the transversality condition for the tensor h. 



, s v^hi,-(n + £/2+ 7/2) -hi,-(n-z/2+ 1/2) 
[Vihi,](n + j72) = J2 ^ I ^-^1 > l_ZJLJ. = 0, (3.21) 

i 

all the terms are defined at the same location n + j/2. 

In Fourier space, the transversality condition is just ki{n)hij{h) = 0, where k(n) is the 
real momentum in Eq. (2.15) and 

h,,(n) = e'^"-("+*~/2+^~/2)h,,(n + i/2 + j/2). (3.22) 

n 

Because the momentum k(n) is real, we can now build the transverse-traceless projec- 
tion in the standard way (1.4)-(1.5), with 

Aij,im(n) = Pu{ii)Pjm{n) - ^P,j (n)Pi„(fi), (3.23) 
P.,(n) = 5.,-^. (3.24) 



hj^{n) = Aijj^ih)him{n) (3.25) 



The projected field 

then satisfies obviously all the requirements: 

1. Symmetry: Aj^.^hi.^ = Aj^- ^^h^^ if = h^^. 

2. Reality: A^;^(n) = Aij-;^(-fi). 

3. Tracelessness: h-u^im = 0. 

4. Transversality: kihij^im = 0. 

5. Maximality: Aij^imi^im = hjj for any symmetric, transverse and traceless hij. 

since the definition of (3.23) is the same as in Section 3.1. 

In actual lattice simulations, it is easier to use a tensor hj^^ defined on lattice sites. 
However, now that we have obtained the projector, we can shift the field to lattice sites by 
a simple translation 

hj;^{n) = hT.T(„ + i/2 + j/2) , (3.26) 

and derive the form of the equivalent projector in that formulation. The Fourier transforms 
of hj^^ and h^""" are related by 

n n 

= e-f ("»+"^)hT.T(fi). (3.27) 

This implies that 

/iT.T(fi) = e-f (^•+^^-)A,,-,„(n)e+t(«^+"™)/,,^(n), (3.28) 
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and therefore, the equivalent projector Aijjm for hjj^ defined on lattice sites is related to 
Ajj,Zm, by 

K,,Un) = e-^(^'+^^)%,„(fi)e+^("'+"'") (3.29) 

Noting the relation between Pim defined as a function of the complex lattice momentum^ 
keff , see Section 3.2, with Pjj defined previously as a function of the real momentum char- 
acteristic of the symmetrized derivative, 



■ IT ~ 



Pi,{n) = e-*iv"»p,^.(n)e+*iv«. , (3.30) 
then we can write the new projector as 

Ai,y^(fi) = Pii{n)Pjmin) - ^e^("'-"^)Pij(fi)Pi™(fi), (3.31) 

P, . % - (3.32, 

At the same time, the shift of coordinates also turns Eqs. (3.19)-(3.21) into a less symmetric, 
but equivalent set of conditions 

Tracelessness : ^ hj^{n — i) = 0, 

i 

Transversality : [/i^'^(n) — hjj^{n — £)] = 0, 

i 

Symmetry : hji^{n) = hj^{n) (3.33) 

Note that Eq. (3.31), together with Eq. (3.32), define a projector that guarantees maxi- 
mality and reality of the Fourier transform of the projected d.o.f. hij{h) = Aij^im{n)him{n), 
as well as transversality, tracelessness and even-symmetry, in the way stated in Eqs. (3.33). 
To be more precise, Eqs. (3.31), (3.32) will guarantee all these conditions with respect 
to the forward/backward derivatives defined in (2.2), which is the one that we will im- 
plement in lattice simulations. For other asymmetric derivatives with lattice momenta 
/ceff,j(n) = e+*^.'Peff,i(n) , (fj G [0,27r), Peffj(n) € Re+, but with e^^'^^ / e+i^"J, the 
projector reads 

A^,Vm(n) = e+^(*'-*')e+^(*^-*-) (^Pu{ii)Pjm{n) - ^e'^^^^-^^^ P,j{n)PUh)^ , (3.34) 

with ^i = ifi — j^hi and Pij{ii) defined as in Eq. (3.32). Of course, in the case of 
the forward derivative, e"*"*"^^ = 6"*"*^^"^^ so $j = V i, and then eq. (3.34) reduces to 
eq. (3.31). Eq. (3.34) will guarantee in general the reality and maximality conditions, the 
even-symmetry and tracelessness defined as in eqs. (3.33), and transversality with respect 
to the lattice derivative defined through the momentum keSj{n) = e"''*'^Jpeff,j(n). 



^Here we refer to the complex momentum fccfr.i of the asymmetric derivative Vi, which differs from the 
real momentum Poff.i of the symmetrized version of Vi, just by a complex phase (pi, i.e. fccft,i/Pcff,j = e"'""''' . 



For the forward derivative this is just e'^^'^' = e 
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4 Comparison of the GW spectra obtained with different projections 



In this Section we will show the time evolution of the the GW spectra during (p)reheating in 
different inflationary models. For each model, we will superimpose the spectra obtained by 
taking Ajjym(n) in Eq. (1.4) evaluated at k = n/ciR (i.e. the continuum-based projector), 
together with the spectra obtained by using the projectors defined in: 1) Eq. (3.2) evaluated 
at the momenta of neutral derivatives, 2) Eq. (3.12) evaluated at the lattice momenta of 
forward derivatives, and 3) Eq. (3.31) evaluated also at the lattice momenta of forward 
derivatives. Thus, we will show explicitly how the several procedures defined to extract 
the GW spectra, compare with each other at different times. We will repeat this exercise 
for different chaotic and hybrid models of inflation. Before describing the results, however, 
we must define the GW spectrum. 

The energy density of a homogeneous and isotropic GW background is described in 
the continuum by [21] 

1 1 



327rG 



/ (ix/iij(x,t)/iij(x,t) 
Jv 

(4.1) 



with (...) a spatial average over a sufficiently large volume V encompassing all the relevant 
wavelengths of the background. In the limit V^^^ — t- oo, the spectrum in the continuum 
(per logarithmic interval) of GW is obtained as 

dpGW f dQk 



I ^h,,ik,k,t)h*^{k,k,t) (4.2) 



dlogk {iTT)^GV J in 

where di^k represent a solid angle element in k-space. 

In the lattice, we must be more careful since clearly we cannot consider the infinite 
volume limit. Assuming the volume of the lattice [V = (Ndx)^] encompasses sufficiently 
well the characteristic wavelengths of the simulated GW background, then it follows 



32ttG iV3 

n 

n n n' 

= ^]^E^r(f^'^)^r(^'0' (4-3) 

n 

where we have used Yin e^*'^™'^^'^"""'-*" = A^^ 5{h — h'). Binning the momentum-lattice in 
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spherical layers of radii |n| and width An, R{n) = {h' / |n| < |n'| < |n| + An}, then 

|n| n'ei?(fi) 

= ;^T^l^V47r|ii|2/ATT(|ii|,t)ATT*(|~| 
|n| 

= Y.{jJi^s ^^(l-l) (/^r(|nNO^r(|nM))^(,)} Alog^ (4.4) 

where A;(n) = km, Alogfc = Ity = ^vr/L, L = Ndx, and (AT.T(|n|,f)/iTT*(|fi|^^)\ 
is an average over configurations with lattice momenta n' G [ |n|, |n| + (5n). 

The spectrum of GW in a lattice of volume V = L^, is therefore defined as 



dpGw \ _ A;^(|n 



dlogkj ^ ' {4tt)^GL3 



(n) = Jri^ ( dx3/iT.T(|fi|,t) dx^/i^:^dn|,t) )^^_^ (4.5) 



„3i,TTn^l 



i?(n) 



In the continuum limit, one identifies'^ DFT{f{n)dx'^} — )• CFT{f{x.)}, and thus, ex- 
pression (4.5) matches the continuum expression (4.2), as it should be. Expression (4.5) 
highlights that the natural momenta in terms of which to express the lattice GW spectrum, 
is the continuum one k = nA^iR, and not any of the lattice-ones defined from the choice of 
a lattice derivative. 

Having defined the appropriate discretized spectrum of GW, we will now show numer- 
ical results from lattice simulations of several scenarios. The key question here will be to 
show the difference in the GW spectra when considering the different projectors defined 
in the previous Section. In particular, we will consider chaotic and hybrid models of infla- 
tion, since the details of the GW production during reheating in those scenarios have been 
studied extensively in the recent years. 

4.1 Transversality check in lattice simulations 

First we define Dififj^ as the sum of hj^ over the lattice sites involved in a particular 
choice of a derivative scheme Vhij. For example, 

Vihf^{w) = y^(/i^'^(n) + hJ^{TL\ + i))/5x, for forward derivatives (4.6) 

i 

VihJ^{n) = '^{hj^{n + + hj^{n - i))/{26x), for neutral derivatives (4.7) 

i 

In the left panel of Fig. 1 we show the evolution in time of the dimensionless ratio 



^One is forced to make this identification since the Parseval theorem in the continuum, i.e. the fact 
that {2Tvf J d^x/2(x) = / d^k !/(k)p, gets translated into {2Tvfdx^ Y^,^ f{x.) = dk^ \dx^f{^)?' in the 
lattice, with dx the lattice spacing and dk = fciR. In the text we refer to DFT and CFT as the discrete 
and continuous Fourier transforms respectively. 
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Figure 1. Left: The time evolution of the degree of transversahty, S{t), obtained for the neutral 
and forward derivatives, both for the continuum- and lattice-based projectors. In the latter case, 
the outcome is clearly limited only by the round-off machine errors, S{t) ^ 0(10^^^), while in the 
former case, S{t) for the continuum can be more than 10 orders of magnitude larger. Right: The time 
evolution of X{t) for the same projectors used in the left figure. Here the degree of transversahty 
is well achieved for all cases, including the continuum-based projector. Note that these plots were 
obtained for a chaotic model, jtf)'^ + ^g'^(px^i with \/g^ = 120. In other models of (p)reheating 
the curves look very similar, with amplitudes of the same order of magnitude. 



for different derivative schemes, where | • • • | is the absolute value and (• • • ) represents an 
average over all the lattice points. Obviously the value chosen for j here is irrelevant. The 
evolution of S{t) gives an idea of how well the transversahty condition VihJ^^ is preserved 
in the lattice. It is clearly appreciated that with the old-continuum projector Ajj /^j the 
transversahty is not that well achieved. Using the projectors which guarantee transversahty 
with respect to the neutral and forward derivatives, Eqs. (3.2), (3.12) and (3.31), we see 
that transversahty is very well preserved (as it should be, by construction), since 5{t) is of 
order 0(10"^^). Such amplitude simply represents the unavoidable round-off errors of the 
machine-precision. Thus, it is very clear that the new lattice-based projectors give raise to 
a well defined notion of transversahty for hj^^ in the lattice. 



We also define the quantity 



m . fy^ (4.) 



and plot it in the right panel of Figure 1, as obtained for all the same projectors used in 
the left panel of the same Figure. As expected, for all cases the degree of tracelessness is 
also very high, only limited again by round-off machine errors. In summary, left and right 
panels of Fig. 1 demonstrate explicitly, and very clearly, that all TT-projectors defined 
in Section 3, effectively filter correctly in the lattice the transverse and traceless d.o.f. of 
two-rank symmetric tensors. 

4.2 GW spectra in the lattice 

Next we will discuss how the new lattice-based projectors modify the GW spectra as 
compared to the spectra obtained with the old continuum-based projector. As we show 
explicitly in Figs. 2, 3 and 4, the spectra of GW in different models is only modified in the 
large-momenta region, i.e. in the ultraviolet (UV) tail. The infrared (IR) features at low 
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momenta, including the shape and amphtude of the spectra, and the position of the peak, 
are not modified by such UV distortion. The UV region corresponds precisely to those 
modes for which the GW spectral amplitude should be exponentially suppressed, if the 
GW spectrum is to be trusted. This is because only the IR modes are excited initially via 
exponential instabilities during (p)reheating [22, 24-26], whereas the UV tail of the spectra 
simply grows by scattering and turbulence [23], see for instance Ref. [3-5] for details. A 
similar behavior occurs also in the context of gauge fields [7, 27-31]. 

That the overall shape and final amplitude^ of the GW spectra does not change much 
when using the new lattice-based projectors, might seem at first sight surprising, given the 
fact that the degree of transversality changes several orders of magnitude when replacing 
the continuum-based projector by the lattice-based ones. However, the lattice-momentum 
/cefr,i(n) from which the lattice-based projectors are made, only differ significantly from the 
momentum used to build up the continuum-based projector, ki = fij/ciR, for the highest 
nj's. For instance, Re{A;^^(n)} = k^g -{h) = sm{2Trhi / N) / dx « /cmnj + 0{2TThi/N)^ as 
long as hi/N < 1/2tt. Thus, as long as nj is not close to the UV boundary of the Fourier- 
lattice hi = zizN/2, and since the Pij operators from which A^j^i^'s are made are quadratic 
in fcefF,!, the difference between the continuum- and the lattice-based projectors can only 
be proportional to the difference |fceff^i(n)p — |/ciRnjp, i.e. 

\A^°U^) - A^5:L(ri)l ~ 0(|keff(ri)P - \hnn\') ~ 0{27r\n\/Nf . (4.10) 

From this point of view, the GW spectra obtained with the lattice-based projectors are 
not expected to differ much from the spectra obtained with the continuum-based projector. 
In particular, in the IR region, say |n| < A^/4, they should be pretty much coincident, the 
better the smaller |n|. Of course, this IR reasoning is still not enough to conclude that 
the GW spectra with continuum- and lattice-based projectors will not be very different. 
As small as it might be such difference, if spurious non-TT modes were incorrectly filtered 
in the continuum-based projector, the difference in amplitude of the two spectra could be 
enhanced during the dynamical evolution of the fields responsible for the GW production. 
That is why implementing in a lattice code the new lattice-based projectors is fundamental 
in order to check whether it makes a difference or not. In Figs. 2, 3 and 4 we quantify this 
aspect, showing the outcome of numerical simulations in which the TT d.o.f. are filtered 
out with the different projectors discussed in Section 3. 

These figures show very nicely the IR aspect just mentioned. Independently of the 
projector, they all tend asymptotically to the same shape in the IR region, the GW spectra 
coincide in shape and amplitude at every step of the evolution, independently of the model 
analyzed. However, the differences in the UV region can be more noticeable, depending on 
the model. But this also depends on the TT-projector used. For instance, in the left panels 
of Figures 2, 3 and 4, we compare the GW amplitude as obtained with the continuum- 
based projector and with the real lattice-based one defined in Eq. (3.2). The Discrepancies 
in all models considered are just a factor 0{1) in the final amplitudes reached, as seen in 
the left panels of all the Figures. Thus, we can conclude that the difference in the UV tails 
are totally irrelevant in this case. 

Let us look now at the right hand side panels of the same Figures, where we compare the 
GW spectral amplitude between the spectra obtained with the continuum-based projector 

*The GW production becomes inefficient in all these models of (p)reheating when the fields enter into 
the turbulent regime, see [3-5] for details, so the spectrum amplitude stops growing and saturates to a 
constant and final shape. 
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Figure 2. The amplitude of the evolution of the GW spectra as obtained with the different 
projectors defined in Section 3, during (p)reheating in the model j4>'^, with A = 10"^'^ and no 
coupling to other fields. Left: Here we compare the spectra obtained with the continuum-based 
projector versus the one obtained with the real lattice-based projector defined in Eq. (3.2). One can 
appreciate that the IR part of both spectra are identical during the whole evolution, whereas the 
UV region shows some difference in amplitude. Such difference smooths out at the end, when the 
different spectra saturate to their final amplitude. Once the GW production becomes inefficient, and 
the spectra of GW does not grow further, the difference among the UV tails between the two spectra 
are simply a factor 0(1). Right: Here we compare the spectra obtained with the continuum-based 
projector versus the ones obtained with the complex lattice-based projectors defined in Eq. (3.12) 
and Eq.(3.31). Again one appreciates that the IR part of both spectra are identical during the whole 
evolution. The difference in the UV region when the spectra saturate to their final amplitude, 
however, is just a factor 0(10), at least at some of the peaks of such UV tail. Note that such 
discrepancy only appears in a region in fc-space where the amplitude of the GW spectrum is already 
a factor 0(10^^) smaller than the maximum amplitude. In any case, the difference between the 
two complex lattice-based projectors is only a factor of order one. 

and with the complex lattice-based projectors defined in eqs. (3.12), (3.31). There we see a 
more noticeable difference in the UV regions. In particular, we observe that in large range 
of momenta, the amplitude of the spectra obtained with the lattice-based projectors can 
be of the order 0(10) bigger than the amplitude for the continuum-based projector. 

In the model jcp'^ with no other couplings to secondary fields^, it is worth noting that 
such discrepancy only appears at a region in fc-space where the amplitude of the GW 
spectrum is already suppressed a factor 0(10"^) compared to the maximum amplitude. 
Therefore, in that respect, it is still a marginal discrepancy. 

In a chaotic model ^(/)^ -|- ^x^'P'^i with (j) the inflaton and x just another field, the 
difference in the UV region is however more visible, see right panel of Figure 3. The 
difference is appreciable at scales in which the GW spectra begin to fall off exponentially, 
but it is not yet suppressed, as compared to the maximum. Such a difference in the 
amplitude of the GW spectra is appreciable by eye, however it only represents an overall 
shift of a factor 0(1) in the location of the scale where the UV tail begins to fall. A similar 
situation arises in the case of a Hybrid model A($^ — f^)^ -|- g^^^x^^ with <I> a waterfall 
field coupled to the infiaton x- In this scenario, see right panel of Figure 4, we find again a 

^In this model, the GW spectra retain the characteristic peaks of the scalar field power spectra, see 
Ref. [4] for more details. 
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Figure 3. The amplitude of the GW spectra obtained with the different projectors defined in 
Section 3. The plots correspond to the evolution of (p)reheating in a Chaotic coupled model 
j(f>'^ + ^-x^^^j with \/g^ — 120. As in figure 2, we can appreciate that the IR parts of the spectra 
are identical during the whole evolution, both in the left and right panels. The discrepancies in the 
UV region are only a factor 0(1) when comparing the output obtained with the continuum-based 
projector versus the one with the lattice-based real projector of Eq. (3.2), see left panel. When 
comparing the GW spectra obtained with the continuum-based versus the complex lattice based 
projectors of eqs. (3.12), (3.31), we see that the difference is more pronounced, a factor 0(10), 
around the scale at which the UV tail begins to fall exponentially. Nevertheless, the difference in 
the UV between the two outputs for the different complex projectors, are still less than a factor 3. 



discrepancy between the lattice-based and the continuum-based projectors at scales where 
the spectra is about to fall off exponentially. Such discrepancy represents nevertheless 
again, only a shift of a factor 0{1) of the scale where the UV tail of the GW spectrum 
begins to fall. 

In any case, both in the hybrid and coupled chaotic models, the difference in the GW 
spectral amplitude between the output obtained with the two complex lattice-based pro- 
jectors (3.12) and (3.31), amounts only to a factor 0{1). So the discrepancies in amplitude 
among spectra obtained with the two complex projectors are again marginal. 

We can conclude that in all the different models considered, the final amplitude, the 
spectral shape of the IR region and the overall shape of the spectra, are not significantly 
modified. There is no leak of scalar modes into the amplitude of the GW spectra even if this 
is obtained with the continuum-based TT-projection, as done repeatedly in the literature. 
The GW spectra show some difference in the UV region when the spectra are extracted 
with the complex lattice-based projectors (3.12) and (3.31). However, those projectors 
present certain caveats, as discussed in Section 3, precisely in the the UV region. When 
comparing the GW spectra obtained with the lattice-based real projector (3.2) and the 
Continuum-based one, the spectral difference in the UV region is indeed really marginal. 
In this latter case, the differences in amplitude in all models considered are just a factor 
0{1), and only show up in the UV region, where the large- momenta tail of the spectra is 
already exponentially suppressed. 
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Figure 4. The amplitude of the GW spectra obtained, for a Hybrid model, with the different 
projectors defined in Section 3. The plots correspond to the evolution after a period of Hybrid 
inflation described by the model A($^ — v'^Y + .g^<i>^x^, with v — 10^^ GeV and A — 2g^. As in 
previous Figs. 2 and 3, we can appreciate that the IR parts of the spectra are again identical during 
the whole evolution, both in the left and right panels. When comparing the output obtained with 
the continuum-based versus lattice-based real projector of Eq. (3.2), see left panel, the discrepancies 
in the UV region are only a factor of order one, so again unimportant. However, when comparing the 
GW spectra obtained with the continuum-based versus the complex lattice-based projectors (3.12), 
(3.31), the difference in the UV is more noticeable. The discrepancy is a factor 0(10) around the 
scale at which the UV tail begins to fall exponentially. Nevertheless, the difference in the UV 
between the two outputs for the different complex projectors, are still less a factor three, similar to 
the Chaotic coupled scenario. 



5 Conclusions and discussion 

In this paper we try to respond to the criticisms made in Ref. [17] with respect to the 
vahdity of the projectors used in numerical simulations of gravitational wave production at 
(p)reheating. It was pointed out that the usual procedure employed in order to obtain the 
Transverse- Traceless part of metric perturbations in lattice simulations was inconsistent 
with the fact that those fields live in the lattice and not in the continuum. It was claimed 
that this could lead to a larger amplitude and the wrong shape for the gravitational wave 
spectra obtained in numerical simulations of (p)reheating, due to the leakage of scalar 
modes into the tensor (GW) modes. In order to address this issue, we have developed a 
consistent prescription in the lattice for extracting the XT part of the metric perturbations. 
We have defined a general complex XT projector based on lattice momenta, as well as a 
real projector in the lattice. All these projectors satisfy the required symmetry properties 
associated with gravitational wave amplitudes. 

We then run specific numerical simulations of GW production at (p)reheating with 
the implementation of the various projectors, and demonstrate explicitly that the GW 
spectra obtained with the old continuum-based XX projection only differ marginally in 
shape with respect to the new lattice-based projectors. Xherefore, we have been able to 
answer the criticisms of Ref. [17] by showing explicitly that the numerical results obtained 
with the lattice-based projectors do not change appreciably with respect to those with the 
continuum-based projector. We thus confirm that all previous results in the literature, 
concerning the spectra of GW coming from lattice simulations of (p) reheating, should be 
trusted to the extent to which such simulations are trusted (i.e. within lattice artifacts' 
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model / projector 


continuum 


lattice real 


latt. complex 


latt. general 


Chaotic Mixed 


7.125 X 10"'^ 


7.15958 X 10-^ 


7.42083 X 10-*^ 


6.46379 X 10-« 


Chaotic Pure 


1.8772 X 10-*^ 


1.87652 X 10"** 


1.88079 X 10-*^ 


1.86935 X 10-« 


Hybrid 


4.6009 X 10-*^ 


4.58099 X IQ-** 


5.15501 X 10"*^ 


5.70762 X 10-« 



Table 1. The total energy density in gravitational waves in units of the total energy density for the 
different models and computed using the various projectors. Note that they differ by less than a 
few percent, except in the hybrid case that can reach 20% between the continuum and the complex 
lattice projectors. 



effects and time evolution). Introducing different lattice momenta in the TT-projector, it 
only gives rise to differences in the (exponentially suppressed) spectral amplitudes in the 
UV, or at most to small shifts (or order unity) in the scale where the spectra begin to fall 
exponentially. The overall shape, frequency and total amplitude do not change. 

Finally, in order to quantify by a single number the deviations induced by the choice of 
projector, to assess the validity of the lattice projector approximation, we have computed 
the fraction of the total energy density in GW to the total energy density available during 
(p)reheating. We have included those ratios in Table 1. Assuming that all the energy at the 
end of inflation went into radiation at reheating, then this fraction represents the fraction 
of GW to radiation today, since they both redshift equally during the subsequent evolution 
of the universe. Therefore the quantity that appears in the table must be multiplied by 
^rad^^ = 3.2 X 10~^ in order to give the observable quantity today, f^GW^^) which could 
eventually be measured in gravitational wave observatories. Note that the entries in Tab. 1 
differ by less than a few percent, except in the hybrid case that can reach 20% between 
the continuum and the complex lattice projectors. These numbers reinforce again the idea 
that the requirement of using the lattice-based projectors does not invalidate the previous 
results found in the literature on GW production in lattice simulations. 
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